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Abstract. The prediction of neutron stars properties is strictly connected to the em¬ 
ployed nuclear interactions. The appearance of hyperons in the inner core of the star is 
strongly dependent on the details of the underlying hypemuclear force. We summarize 
our recent quantum Monte Carlo results on the development of realistic two- and three- 
body hyperon-nucleon interactions based on the available experimental data for light- and 
medium-heavy hypernuclei. 


1 Introduction 

Neutron stars are among the most compact and dense objects in the Universe, with typical masses 
up to ~ 2 Mq and radii R ~ 10 km. Their central densities can be several times larger than nuclear sat¬ 
uration density, po - 0.16 fm“^. At densities larger than ~ 2po, when the nucleon chemical potential 
is large enough, the conversion of nucleons into hyperons might become energetically favorable, and 
Pauli blocking would prevent hyperons from decaying by limiting the phase space available to nucle¬ 
ons. This would lead to a reduction of the Fermi pressure exerted by the baryons and to a softening 
of the equation of state (EOS). As a consequence, the maximum mass determined by the equilibrium 
condition between gravitational and nuclear forces would be reduced. 

Currently there is no general agreement (even qualitative) among the predicted results for the EOS 
of strange baryonic matter and the maximum mass of a neutron star (NS), that violates the constraints 
given by the recently measured value of ~ 2 Mq of the millisecond pulsars PSR J1614-2230 [1] 
and PSR J0348H-0432 [2], see Ref [11]. This apparent inconsistency between NS mass observations 
and theoretical calculations is a long standing problem known as the hyperon puzzle. It has to be 
ascribed to a combination of an incomplete knowledge of the forces governing the system and to 
the concurrent use of approximate theoretical many-body techniques. Its solution requires a more 
accurate theoretical analysis and a thorough experimental investigation of the hyperon-nucleon force 
in a variety of systems, ranging from light to medium and heavy hypemuclei. 

We employ a quantum Monte Carlo computational technique to tackle the general problem of 
the hyperon-nucleon interaction. The algorithm, auxiliary field diffusion Monte Carlo (AEDMC), 
has been introduced by Schmidt and Eantoni [3] as an extension of the usual diffusion Monte Carlo 
method to deal in an efficient way with spin/isospin-dependent Hamiltonians. It has been successfully 
applied to solve the many-body problem in a non perturbative fashion for strongly correlated systems, 
both in the standard [4-7] and strange nuclear sector [8-11]. 
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2 Hypernuclei: from few- to many-body systems 

As reported in Refs. [8-10], an accurate analysis of the A separation energy B\ of light- and 
medium-heavy hypernuclei has been recently carried out using a phenomenological hypernuclear 
interaction analog to the Argonne-Illinois nucleon-nucleon force [12-19]. The AN two-body compo¬ 
nent of this interaction has been fit to the existing Ap scattering data. The ANN three-body sector 
instead has been re-fit to the available A separation energies of single-A hypernuclei over a wide mass 
range by means of AFDMC calculations. It has been shown that the repulsive nature of the three-body 
hyperon-nucleon interaction is the key to satisfactorily reproduce the ground state properties of hyper¬ 
nuclei from light- to medium-heavy within a unique theoretical framework. Fig. 1 shows the AFDMC 
results for when using the complete AN + ANN force. The agreement with experimental data is 
good over all the mass range investigated and for the A particle in different single particle states. For 
more details see Refs. [8-10]. 

In the light sector, a two-body AN charge symmetry breaking interaction has been also em¬ 
ployed [9, 13]. As shown in Tab. 1, AFDMC results reproduce the hyperon separation energies and 
energy splitting of the mirror hypernuclei and ^He, according to the old emulsion data [20]. 
is weakly bound compared to deuterium and the neutral hypernucleus ^n is predicted to be unbound, 
in agreement with other theoretical results [21-23], contradicting the possible experimental evidence 
reported in Ref. [24] . 

In order to validate the phenomenological approach adopted in quantum Monte Carlo calcula¬ 
tions for hypemuclei, a benchmark among different few-body methods has been recently proposed. 
In tab. 2 we report preliminary results for ^He for two different exact methods, AFDMC and Non- 



A-2/3 

Figure 1. Measured (see Ref. [10] for the complete list of references) and computed A separation energies as a 
function of . Results for the A particle in different single particle states are shown. 
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Table 1. AFDMC and experimental [20] A separation energies for the lightest hypemuclei. In the case of A = 4 
the separation energy difference AB^ is also shown. Evidence for the neutral hypernucleus ^n has been reported 
in Ref. [24], but no Ba value has been reported. All energies are in MeV. 


System 

Ba 

B7 

ABa 


Hn 

unbound 

unknown 




0.23(9) 

0.13(5) 



4 He (0") 

1.95(9) 

2.37(9) 

2.04(4) 

2.39(3) 

0.42(11) 

0.35(6) 


symmetrized Hyperspherical Harmonics (HSHH) [25]. In both strange and non-strange sectors the 
employed potentials are at two-body level only. Results agree among different methods whether using 
a soft (Minnesota [26]) or hard (Argonne V4' [27]) NN potential. 

Table 2. Binding energies and A separation energies for ^He for different combinations of nucleon-nucleon and 
hyperon-nucleon potentials and different methods. All energies are in MeV. 


^ He 

AFDMC 

NSHH 

NN Minn [26] + AN 

E -37.69(8) 

E 

-37.77(10) 


Ba 6.95(9) 

Ba 

6.99(10) 

NN AV4' [27] -H AN 

E -39.46(12) 

E 

-39.54(10) 


Ba 6.70(16) 

Ba 

6.84(10) 


3 Neutron stars: from hypernuclei to the infinite medium 


Precise AFDMC calculations have been performed to first assess the parameters of the interaction 
from measured hyperon separation energies, and then to extrapolate the results to the case of neutron 
matter with strangeness. However, parametrizations of the potential predicting relatively small differ¬ 
ences in the A separation energies of hypernuclei give very different results for the properties of the 
infinite medium [11]. The resulting EOS spans the whole regime extending from the appearance of a 
substantial fraction of hyperons at ~ 2po to the absence of A particles in the entire density range of the 
star. This yields a sizable effect on the predicted NS structure. Therefore, the derivation of realistic 
hypernuclear potential models is of paramount importance to properly assess the role of hyperons in 
NSs and reconcile theoretical predictions with astrophysical observations. This demands a precise 
and systematic experimental investigation of hypernuclear properties over a wide range of masses. 
For instance, a recent study of the isospin dependence of the present three-body hyperon-nucleon 
force shows the difficulty in extracting the information on the Hamiltonian from currently available 
experimental information on hypernuclei [10]. 
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